Ab initio investigation of structural and magnetic properties of Fe|ZnSe(001) 
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We investigate 1, 2, 3, and 4 monolayers of iron grown on the (001) surface of ZnSe using spin 
density functional theory, plane-wave expansion, and ultrasoft pseudopotentials. An ideal structure 
of the interface is plausible due to the very small lattice mismatch between both materials. However, 
we observe dramatic deviations from this idealized configuration. The most distinctive variations 
arise in the case of a single monolayer. Here a terminating Zn layer is buoying upwards. The thicker 
the Fe layer the weaker are the relaxations in the iron region. Nevertheless, in the case of the Se- 
terminated surface, the structure changes significantly inside the semiconductor near the interface. 
The profiles of the electronic density and the magnetization perpendicular to the surface are quite 
sensitive to the change from the ideal to the real structure. 

PACS numbers: 68.35.-p, 75.50.Bb, 75.70.Ak, 81.05.Dz 
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I. INTRODUCTION 

Magnetoelectronicsi and spintronicsi are rapidly grow- 
ing fields in applied solid state physics. Since the 
proposal of a spin field-effect transistor by Datta and 
Das£ ferromagnet-semiconductor (FM-SC) heterostruc- 
tures have become very interesting candidates for the 
realization of spin dependent injection of electrons in 
a semiconductive material. Therefore, theoretical^^ 
and experimenta&S*^ effort has been made. Towards 
this goal, calculations were made for FM-SC-FM het- 
erostructures, while thin Fe films epitaxially grown on 
semiconductor surfaces were investigated experimentally. 

The first Fe|GaAs(001) systems- had shown magnetic 
dead layers. This phenomenon of Fe layers at the in- 
terface between the two materials with nearly vanishing 
magnetic moment can be prevented by lowering the tem- 
perature of the substrate to room temperature during 
the film growth. 10 Higher temperature leads to a mixing 
of the materials at the interface which sensitively influ- 
ences the magnetic properties. Thus, structure plays an 
important role. 

Nevertheless, the theoretical investigations of FM-SC 
heterostructures usually assume an ideal interface struc- 
ture. This is suggested by the very small mismatch of the 
lattice constants between diamond-structure elemental or 
zincblende-structure compound semiconductors and of a- 
Fe. Therefore, four Fe cubes excellently fit on a single 
fee cube of the semiconductor material. The advantage 
of these computations is the possibility of considering 
the contribution of orbital magnetic moments. This is 
done within the LKKR-methodii which provides the op- 
portunity to take into account the relativistic spin-orbit 
coupling. However, the possibility of sensitive reactions 
of the total magnetization to minor structural changes is 
fully neglected. 

In this article we show that these idealizing assump- 
tions are not realistic in many cases. Therefore, we deter- 
mine the equilibrium structure of the observed systems 
and compare the ideal and real density profiles of the va- 
lence electrons and of the magnetization. We will see that 



the deviations from the ideal behavior transcend small 
quantitative corrections even for the structural results. 
Because the orbital contribution can be shown^ to con- 
tribute less than about 5% in the case of the ideal struc- 
ture, we consider only the spin magnetization as the main 
contribution to the total magnetization. This quantity 
is directly available from spin density functional theory 
and thus easier to obtain. Furthermore, the small con- 
tribution of the orbital moment can be neglected, if we 
expect a sensitivity of the magnetic properties to struc- 
tural changes, which should affect the magnetization pro- 
files by more than a few percent. 

We choose Fe|ZnSe(001) as an example for the whole 
class of possible systems, because it is of experimen- 
tal interest^ and the pristine semiconductor surface 
shows much smaller reconstructions than any III-V 
semiconductor ii2ii^ A II- VI material as ZnSe can satisfy 
the well established but not rigorous counting rulesi^i^ 
by building simple dimers at the (001) surface. Hence, 
the often used GaAs(OOl) surface is much more compli- 
cated. 



II. METHOD 

Density functional theory (DFT)i^ii& provides a reli- 
able ab initio method for computing the ground-state 
energy and the associated one-particle density of an 
electronic many-particle system with given external po- 
tential. Thus, it is possible to determine the ground- 
state configuration of ions in a chemical system, such as 
molecules or crystals. Due to the ferromagnetic charac- 
ter of iron we need to consider the spin degree of free- 
dom. This possibility is offered by spin density functional 
theory (SDFT)i 17 i 18 i 19 Within this generalization of the 
original DFT the ground-state energy is a functional of 
the density as well as of the magnetization. 

In this work, we use the ab initio package VASR2i2!i2£ 
for the numerical implementation of DFT and SDFT. 
Due to the three-dimensional lattice symmetry of the 
bulk systems the wave functions of the valence electrons 
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TABLE I: Cubic lattice constants of iron and zincselenide as results of LDA an d GGA calculations, bulk moduli A'o and 
the according pressure derivatives K' . The experimental data is taken from Refs. I20j2lll22l The lattice constant for 7-Fe is 
measured at a temperature of 916 °C, the other data corresponds to room temperature. 



are expanded in in terms of plane waves. The core 
states are accounted for by ultrasoft pseudopotentialsSi 
similar to Vanderbilt's potentials ^1 which are no longer 
norm conserving. The exchange and correlation energy 
is treated within the local-density approximation (LDA) 
and the generalized-gradient approximation (GGA). In 
all LDA functionals the parameterization of Perdew and 
Zunger 2 >2. to Monte-Carlo data of Ceperly and Alder— is 
used; in the GGA functionals the proposal of Perdew and 
Wang (1991) 3 £ is adopted. 

Using these techniques we are able to determine the 
minimal-energy lattice constants of zincblende ZnSe, 
body-centered cubic (bcc) a-Fe, and face-centered cubic 
(fee) 7-Fe. Therefore, we compute the energy per unit 
cell as a function of the cubic lattice constant. The cut- 
off energies for the expansion in in terms of plane waves 
are 209. 5eV for ZnSe and 297eV for Fe. The momentum- 
space integrations are done by summing over a grid of 
special points according to Monkhorst and Pack— Here 
we use a 10 x 10 x 10 mesh for ZnSe with two atoms per 
unit cell and 20 x 20 x 20 meshes for both modifications 
of iron each with one atom per cell. The results of the 
Murnaghan— fits to the results are given in Table [I] The 
lattice constants agree very well with the experimental 
data and with the result of former calculations— iSiSiSi 
Also the bulk moduli and their pressure derivatives, 
which result from the Murnaghan fits, show good agree- 
ment with experiment and computation. Additionally, 
iron is ferromagnetic in the bcc structure. The sponta- 
neous magnetic moment per atom is /j, = 2.04 /zb using 
LDA and /1 = 2.35 [Ab in a GGA computation. Thus, 
the values \i = 2.12 fiQ using LDA and /1 = 2.35 /zb us- 
ing GGA of a former calculation^ can be reproduced. 
These results soundly fit to the experimental^ value of 
/i = 2.23 /ib- However, using LDA 7-Fe turns out to 
be stable in fee and not the actual bcc structure. This 
inconsistence with the experiment is an artifact of the 
local-density approximation— which can be repaired us- 
ing generalized-gradient corrections. 

Surfaces are treated within the periodic-slab method— 
The program package VASP needs three-dimensional lat- 
tice symmetry, which is broken in the case of a surface. 
But it is possible to restore periodicity perpendicular to 
the surface by periodically stacked slabs of material. Here 
we take into account nine atomic layers of the semicon- 
ductor ZnSe covered by one, two, three, or four layers of 



iron. The dangling bonds at the uncovered side of ZnSe 
are saturated with hydrogen atoms. The distance be- 
tween the slabs is about 8.5 A. Tests showed that these 
parameters are a good compromise between computa- 
tional effort and accuracy. More details about the ge- 
ometry of the used supercell can be taken from Sec. IIIII 
and Fig. ^ 111 these computations we use a 6 x 6 x 1 
mesh of special points and an energy cut-off of 237.5 eV. 
The ionic relaxation to the equilibrium structure is de- 
termined using a conjugate-gradient algorithm- 
Ill. THE IDEAL INTERFACE 

The lattice mismatch 7 of a-Fe and ZnSe can be de- 
fined as follows: 

aznSo — 2ap c 

7= W -7, 7- W 

2 (aZnSc + la-pe) 

Comparison with Table ^ shows that this mismatch is 
7 = 1.0% using LDA, 7 = 0.2% using GGA, and 
7 = —1.1% for the experimental values. Thus, the lat- 
tice mismatch is very small, and the idea of an ideal 
structure of the interface is suggestive. For many sys- 
tems of iron grown on the (001) surfaces of zincblende 
semiconductors the mismatch of the lattice constants is 
quite small (experimental values*-* e.g., acaAs = 5.653A, 
oaiAs = 5. 660 A). Hence, many calculation o 4 ^^ are based 
on this idealization. 

Fig. Ogives an intuitive picture of the ideal interface. 
Regarding the pristine surface of ZnSe first, it is evident 
that (001) surfaces always are terminated by only one 
sort of atom, thus we have to distinguish between the 
Se- and Zn-terminated case. Due to the tiny lattice mis- 
match four iron cubes closely fit on a single ZnSe cube. 
Substitution of the lowermost atomic layer of iron by the 
first semiconductor layer results in the ideal structure 
shown in the middle of Fig. ^ In this case the semicon- 
ductor is Zn-terminated. 

But the presented cell is not elementary in a two- 
dimensional sense parallel to the surface. However, an 
elementary cell can be obtained by rotating the structure 
around an axis perpendicular to the surface by an angle of 
45° and cutting off the corners. This is visualized on the 
right side of Fig. ^ by showing the elementary cell used 
in the calculations for the Zn-terminated semiconductor 
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covered by two monolayers of iron. As already mentioned 
in Sec. [HI we choose 9 atomic layers of ZnSe and 1, 2, 3, 
or 4 monolayers of iron. The separating vacuum slab of 
a thickness of about 8.5 A is not displayed. A monolayer 
in ZnSe contains either one Zn or one Se atom in the 
cell, a single monolayer of Fe contains two inequivalent 
atoms. The dangling bonds at the bottom side are satu- 
rated with hydrogen atoms. The positions of the H atoms 
and the lowermost atoms of ZnSe are fixed in the calcu- 
lations; thus, the bond distances between H and Se or Zn 
are constant with <inse = 1-478 A and dnzn = 1-612 A. 
These values are the equilibrium bond distances of the 
corresponding isolated dimers computed with VASP. 

Starting from the ideal structure it is possible to com- 
pute the relaxation of the atoms in the elementary cell. 
If restricted to the direction parallel to the surface nor- 
mal the displacements turn out to be small, and thus 
the structure of the surface remains nearly ideal. How- 
ever, this constraint may lead to possibly unstable equi- 
librium positions. In order to allow relaxation into recon- 
structed equilibrium positions we start from randomly 
displaced positions of the peripheral Fe atoms. The re- 
sults of these computations may differ dramatically from 
the ideal structure. This is reported in the following sec- 
tion. 



IV. RESULTS 

In this section we present the results obtained for Se- 
and Zn-terminated (001) surfaces of ZnSe covered by 1, 
2, 3, and 4 monolayers of iron. In all cases the energy- 
minimizing structure is determined first. For an esti- 
mation of the dependence of magnetic properties on the 
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FIG. 1: Structure of the ideal interface. Here the semi- 
conductor surface is Zn-terminated. 4 Fe cubes fit on a single 
ZnSe cube as displayed on the left. Thus, the structure shown 
in the middle is plausible. The elementary cell used for the 
calculation is presented on the right. 



structure we show the profiles of the spin magnetization 
perpendicular to the interface. The term profile means 
that the corresponding magnitude is averaged over planes 
parallel to the surface. In bulk iron the orbital magne- 
tization is suppressed because of quenching by electric 
crystal fields4i At the surface or interface the crystal 
symmetry is broken; thus, the quenching does not suffice 
any more to disable completely the orbital contribution 
to the magnetization. Nevertheless, in calculations^ con- 
sidering the orbital contribution this appears to be not 
more than a 6% effect. Thus, the mere spin polariza- 
tion is a good indicator for the behavior of the magnetic 
properties when structural deviations from the ideal case 
have to be expected. 



A. A single monolayer of iron 

The most interesting results are obtained for only one 
Fe monolayer on ZnSe(OOl). In particular, the relaxed 
structure deviates tremendously from the ideal case. 
Fig. [21 shows the peripheric region of the covered Se- and 
Zn-terminated surface. In the Se-terminated case the two 
inequivalent Fe atoms obviously have very different ver- 
tical positions. One Fe atom is lowered, which enables a 
bond to a Zn atom of the second atomic layer of ZnSe. 
The other Fe atom is bonded to the uppermost Se and 




Zn-terminated 



FIG. 2: Structure of the Se-terminated (upper part) and Zn- 
terminated (lower part) ZnSe(OOl) surface covered by a single 
monolayer of iron. Only the crucial region near the interface 
with interesting distances and angles is shown. 
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the other Fe atom. Fe and Se together build up a surface 
slab that is coupled to the ZnSe substrate only by an Fe- 
Zn bond. So there is no bond between Se and Zn at the 
two uppermost layers any more. In this slab, the the Fe 
and Se ions form zigzag chains from the front backwards 
(in Fig. 0) with bond lengths di = 2.30 A, d 2 = 2.42 A, 
and Se-Fe-Se bond angle a = 119°. Similar chains run 
from the left to the right with the values d-j = 2.34 A, 
d A = 2.33 A, and = 121°. 

The Zn-terminated surface provides the most amaz- 
ing results. The uppermost Zn layer is buoying upwards 
which can be seen at the lower part of Fig. [21 We obtain 
two layers in the perpheric region. The uppermost layer 
consists of one Fe and one Zn atom in the elementary 
cell with nearly identical heights. The second layer is 
made up of one Fe and one Se atom per cell. Similar to 
the Se-terminated case these atoms build zigzag chains 
sidewards. There the bond distances are ds = 2.44 A, 
d e = 2.42 A, and the Se-Fe-Se angle is = 121°. Again, 
the uppermost Se atom is only bonded to iron, and the 
hole peripheric layer is coupled to the semiconductor by 
an Fe-Zn bond. Both structural results can be inter- 
preted in the way that non-metallic Se forms tight bonds 
to Fe. This affinity is so distinctive that iron can even 
displace the originally terminating Zn atoms towards the 
outermost region. 

An indicator for the stability of the surface structure 
is the energy gain due to the relaxation in comparison 
with the ideal structure. Also the results for DFT and 
SDFT have to be differentiated. Because the relaxation 
was done initially using usual DFT the result of which 
was refined with SDFT afterwards, we have to distin- 
guish between energies from DFT and SDFT. For the 
Se-terminated surface the energy gain is AE± = 0.638eV 
from the small relaxation of the peripheric atom perpen- 
dicular to the interface. Complete DFT relaxation causes 
an increase to AEdft = 1.663 eV by more than 1 eV. 
This value is much smaller using SDFT even after further 
relaxations. The smaller energy gain A^sdft = 0.346eV 
can be interpreted in the following way: DFT only can 
give the energy differences due to structural changes, 
whereas using SDFT this is corrected by the energy due 
to the ordering of the spins. Thus, the ferromagnetic or- 
der in the ideal structure is energetically advantageous 
in comparison with the relaxed structure. Nevertheless, 
the structural contribution to AE outbalances the ferro- 
magnetic one; thus, the energy gain is still much greater 
than the thermal energy k^T ss 25 meV corresponding 
to room temperature. 

In the case of the Zn-terminated surface the respective 
values are AE±_ = 0.125 eV, AE DFT = 2.767 eV, and 
AE SBFT = 0.577 eV. AE± is smaller than in the Se- 
terminated case, because the reduction of the distance 
between the metal Fe and the nonmetal Se yields a higher 
energy gain than a variation of the Fe-Zn bond length. In 
the DFT calculation, AE'dft is much larger than in the 
Se-terminated case. So the energetic advantage of the Zn 
atoms buoying upwards is enormous. On the other hand, 




FIG. 3: Structure and density profiles for one monolayer of 
Fe (2 atoms per cell) on Se-terminated ZnSe(OOl). The ver- 
tical dotted lines mark the positions of the atomic layers in 
the ideal case. At the top the elementary cell of the ideal 
interface is presented. At the bottom the atoms are shown 
after relaxation in two different perspectives. The uppermost 
plot shows the density of the valence electrons, in the mid- 
dle the magnetization is displayed, and in the lowermost plot 
the magnetization is shown on an expanded scale. The solid 
(dashed) lines belong to the relaxed (ideal) structure. 



this dramatic structural variation reduces the energy gain 
due to ferromagnetism more than for the Se-terminated 
surface. But still A.Esdft is much greater than typical 
thermal energies at room temperature. 

The density profiles measure the dependence of elec- 
tronic and magnetic properties on structural variations. 
The profiles for the density of the valence electrons, the 
magnetization, and the magnetization on an expanded 
scale are shown in Figs. |3] and 0] Additionally, the atomic 
positions in the elementary cells are presented for the 
ideal and relaxed structure. The valence-charge density 
has minima for Se layers and maxima for Zn and Fe lay- 
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FIG. 4: Same as Fig. El but for one monolayer of Fe on 
Zn-terminated ZnSe(OOl). 



ers. This is because Se is treated with a pseudopotential 
for 6 valence electrons (4s 2 4p 4 ), whereas for Zn 12 elec- 
trons (3d 10 4s 2 ) are taken into account, although the total 
number of electrons is larger for Se [Z = 34) than that 
of Zn (Z = 30). All the changes of the electron den- 
sity can be interpreted as a shift of the density along 
the direction of the corresponding relaxations. Also the 
magnetization seems to follow the ionic relaxations, es- 
pecially for the Fe layers, although the total magnetic 
moment per cell decreases from 6.32 /xb for the ideal in- 
terface to 4.94 /ie after relaxation in the Se-terminated 
case and from 6.17/iB to 5.34 /iB, if the ZnSe surface was 
Zn-terminated. The exact run of the curves in Figs. El 
and 0] according to the relaxed structure differs from the 
ideal case, which is clear because of the splitting of the 
Fe plane into two layers. 

The lowermost graphics in both, Fig. El and EJ show 
drastic differences between the magnetization before and 
after relaxation on an expanded scale. In the case of the 
Se-terminated surface (Fig. E| the magnetization inside 



the semiconductor even has the direction opposite to the 
total moment. This antiferromagnetic ordering is obvi- 
ously an effect of the structural changes, because for the 
ideal interface the magnetization has the same direction 
inside the semiconductive and the ferromagnetic region. 
An appropriate interpretation is given by the Heisenberg 
models There, the coupling constants Jij are given by 
the overlap of the wave functions of neighboring atoms. 
Because the comparatively strongly localized d-electrons 
dominate the ferromagnetism of Fe, a small change of 
bond lengths or layer distances can invert the sign of 
the coupling constants, which means a transition from 
ferro- to antiferromagnetism. Additionally, the relax- 
ation results in an Fe-Zn bond that does not exist in 
the ideal case. In contrast to Fe-Se bonds, these bonds 
lead to a flip of the magnetization, as can be seen for 
the Zn-terminated interface. For the Zn-terminated sur- 
face (Fig. EJ) a more complicated run of the curves inside 
ZnSe is obtained. Fe-Zn bonds at the interface instead 
of Fe-Se bonds for Se termination cause antiferromag- 
netic effects, even for the ideal structure. However, the 
magnetization deep inside the semiconductor has flipped 
its direction after relaxation, too. The tiny peaks at the 
position of the hydrogen atoms indicate too thin slabs 
or vacuum layers. But incrementing those two parame- 
ters would require a drastic increase of computing time. 
Because we just want to estimate the dependence of mag- 
netic properties on the structure, the chosen accuracy is 
sufficient. 



B. Two, three, and four monolayers 

Now we turn to the case of 2, 3, or 4 monolayers of Fe 
grown on ZnSe(OOl). The structure and the associated 
density profiles are presented in Figs. 03 El and [7| for the 
Se-terminated surfaces and in Figs. H3 El and ITU1 for the 
Zn-terminated case. Here, we omit the presentation of 
the structure like the one in Fig. El in a separate plot. In 
the results for a given termination clear similarities arise, 
thus a common discussion is sensible. 



1. Se termination 

The structures of all three systems are very similar. 
Those are presented beneath the plots of the profiles in 
Figs. El (2 monolayers), El (3 monolayers), and [7] (4 mono- 
layers). The iron slab has nearly perfect bec structure as 
known from bulk a-Fe. Relative relaxations parallel and 
perpendicular to the interface are negligibly small. How- 
ever, the whole Fe region is relaxing towards the semi- 
conductor substrate by several 0.1 A, which is an effect 
of structural changes in ZnSe. The first atomic layer of 
ZnSe does not solely exist of Se any more. The Zn atoms 
of the originally second atomic layer together with the 
terminating Se atoms build a ZnSe layer at the interface. 
So it is possible that the uppermost Se atoms are not 
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bonded to Zn but to Fe, thus the whole surface slab is 
coupled to the substrate by Fe-Zn bonds followed by a 
Zn-Se bond. The distance of the second Se layer, which 
can be regarded as the uppermost layer of the substrate, 
to the first layer of the peripheric part of the slab (the 
ZnSe layer) is clearly greater than the layer distances in 
the ideal structure. So the peripheric part of the slab 
should be coupled weakly to the substrate, and it can be 
displaced easily in the horizontal direction. This would 
implicate a very soft phonon mode. Measuring such a 
mode could confirm our results for the structure. Deeper 
in the semiconductor the structure becomes more and 
more similar to that of zincblende. Beyond the fourth 
atomic layer (second Zn layer) the structure can be re- 
garded as ideal. The energy gain per cell decreases with 
increasing number of Fe layers from 0.451 eV to 0.392 eV 
and 0.279 eV due to the increasingly ideal character of 
the ionic configuration. 

The almost ideal structure of the iron is reflected in the 
nearly ideal profiles of the valence electron density and 
the magnetization in the Fe region, too. The peaks cor- 





FIG. 5: Same as Fig. [3] but for two monolayers of Fe on 
Se-terminated ZnSe(OOl). 



FIG. 6: Same as Fig. |2l but for three monolayers of Fe on 
Se-terminated ZnSe(OOl). 



responding to iron are simply shifted along the direction 
of the ionic relaxation. Only the peak nearest to the in- 
terface is slightly decreased for the case of two and four 
monolayers of Fe. So the total magnetic moment per 
cell reduces from 11.14 /xb for the ideal structure with 
two Fe layers to 10.68 /xb after relaxation, from 16.79 /xb 
to 16.35 /ib for three monolayers and from 21.71 /j,b to 
20.49 /is for four monolayers. The ZnSe layer at the 
interface entails an additional peak in the electronic den- 
sity, because this layer is packed very densely. Towards 
ZnSe this peak is followed by a dip, which clearly shows 
the already assumed weakness of the expanded Zn-Se 
bond. Inside the semiconductor the density corresponds 
to a shift along the relaxation again. The magnetization 
inside ZnSe declines very fast in all three cases. However, 
the direction of the weak magnetization is flipped, simi- 
lar to the case of a single monolayer of Fe (compare Fig. 
|3J). Again different layer distances and changing chemical 
environments result in antiferromagnetic ordering. In all 
three cases the magnitude of the magnetization and the 
penetration depth slightly increase. 
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FIG. 7: Same as Fig. |31 but for four monolayers Fe on Se- 
terminated ZnSe(OOl). 



2. Zn termination 

For the Zn-terminated surfaces a different picture of 
the structure is obtained. Here the deviations from ideal 
behavior preferentially appear in the Fe slab, whereas 
ZnSe remains in the almost perfect zincblende structure. 
Inspection of Figs. 00 E3 and 1101 demonstrates that the re- 
laxations in the semiconductor region are nearly vanish- 
ing. There is no drastic structural change at the interface 
like the formation of a ZnSe layer in the Se-terminated 
case. However, relaxations in the Fe region are obvious. 
Here the Fe atoms are displaced in the direction paral- 
lel to the interface; thus, the bond lengths change, but 
the layer distances do not. For only two monolayers of 
Fe these relaxations are most evident. Actually, the Fe 
atoms at the interface are able to form additional bonds 
within the layer. The deviations from the ideal struc- 
ture seem to become weaker with an increasing number 
of monolayers. In contrast to the Se-terminated interface 
the Fe slab slightly relaxes outwards. 

Because of the smaller extent of the structural mod- 
ification in comparison with the Se-terminated case we 



expect smaller energy gain due to relaxation. But sur- 
prisingly the values 1.027 eV, 0.814 eV, and 0.610 eV are 
much greater than for Se termination. Changing the ionic 
configuration in the Fe slab seems to reduce the energy 
by a greater amount than the relaxations in the semicon- 
ductor do. However, the energy gain decreases with the 
number of Fe layers, which again shows the increasingly 
ideal character of the system. 

The density profiles of valence electrons do not ex- 
hibit any astonishing effects. Because there are no closer 
packed layers or strongly increasing layer distances, only 
the outwards relaxation of the whole Fe region can be 
identified in a shift of the Fe peaks. However, the mag- 
netization even within the iron region extremely deviates 
from the ideal behavior. For the ideal structure the mag- 
netization due to the Fe layer at the interface is strongly 
suppressed in the cases of 2 (Fig. [HJ) and 3 (Fig.^J mono- 
layers of Fe. This result could be interpreted as magnetic 
dead layers claimed in older work. 7 However, the effect 
of nearly vanishing magnetization at the interface turns 




FIG. 8: Same as Fig. [3] but for two monolayers of Fe on 
Zn-terminated ZnSe(OOl). 




FIG. 9: Same as Fig. [3] but for three monolayers of Fe on 
Zn-terminated ZnSe(OOl). 



FIG. 10: Same as Fig. but for four monolayers of Fe on 
Zn-terminated ZnSe(OOl). 



out to be unrealistic: In the case of four monolayers of 
Fe all magnetic moments are fully developed even for the 
ideal structure, and in the case of two and three mono- 
layers the dead layers disappear after the relaxation. So 
we find in agreement with experiment^ and other theo- 
retical work 6 that there are no magnetic dead layers in all 
realistic systems considered here. Due to the difficulties 
in the treatment of the ideal interfaces the total mag- 
netic moment per cell increases for two monolayers from 
7.83 /ib before relaxation to 11.09 /ib after and for three 
monolayers from 10.26 /zb to 16.64 /iB- The known effect 
of a slight decrease of the moment as seen for Se termi- 
nation is obtained again for four monolayers of Fe. Here 
the moment reduces from 21.62/xb to 21.52/xb by the tiny 
value of 0.1 fiB- The difference is obviously smaller than 
for Se termination. Amazingly, the relaxations of the 
Fe atoms parallel to the surface can stabilize the mag- 
netic moment, whereas the decrease of the moment is 
even less pronounced for Se-terminated surfaces, where 
the iron structure remains almost ideal. 

Inside the semiconductor the magnetizations for the 
ideal and the relaxed structure strongly deviate from each 



other in the case of two and three monolayers of Fe. This 
is no surprise because of the dead layers in the ideal struc- 
ture. So the differences are difficult to interpret, the more 
so as the lowermost profile in Fig. 1101 for four monolay- 
ers shows only little difference between the two curves 
corresponding to the ideal and relaxed structure. For all 
Zn-terminated systems including a single monolayer of 
Fe a tendency towards oscillations of the magnetization 
is evident. This is in contrast to Se termination. 



V. SUMMARY AND CONCLUSIONS 

In this work we have presented structural, electronic, 
and magnetic results from spin density functional the- 
ory calculations for 1, 2, 3, and 4 monolayers of Fe on 
Se- as well as Zn-terminated ZnSe(OOl). In spite of the 
tiny mismatch of the lattice constants of a-Fe and cu- 
bic ZnSe we find dramatic structural deviations from the 
ideal interface. Especially for a single monolayer of iron 
the structure is changed radically. But also for more than 
one monolayer we obtain amazing results, so we find a 
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separated ZnSe layer with Zn and Se at Fe positions at 
the interface in the case of Se termination. Otherwise, for 
Zn termination the structural modifications take place in 
the Fe region. 

Summarizing, we have to point out that the structures 
become increasingly ideal with an increasing number of 
Fe monolayers. We see the tendency of tight Fe-Se bonds; 
in every Se-terminated system the outermost Se atoms 
are bonded primarily to Fe. The raised Zn atoms are 
bonded to Fe and form a single weak bond to the sub- 
strate. In the Zn-terminated structures, too, the whole 
surface layer is coupled to the substrate by an Fc-Zn 
bond. 

We find a clear dependence of the profiles of electronic 
density and magnetization on the relaxations. The elec- 
tronic density and the magnetization due to the Fe atoms 
can be explained by a shift of the curves along the direc- 
tion of the relaxation. However, the magnetization inside 
ZnSe often changes the direction in comparison with the 
ideal structure or the total magnetic moment. Although 
the magnitude of this magnetization is small, the flip of 



the sign might affect electronic and magnetic properties 
near the interface. For example a reliable description 
of the tunneling magneto resistance (TMR) requires an 
accurate electronic density of states at the interface 
However, this magnitude might be sensitive to small re- 
laxations due to changing bond conditions. 

Thus, the effect of relaxations should not be ne- 
glected in more realistic calculations of electronic and 
magnetic properties of ferromagnet-semiconductor het- 
erostructures. The usage of the ideal interface in such 
computations can foil ambitious efforts using complicated 
methods. 
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